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ABSTRACT 

We present an analytic expression for the axis ratio distribution of triax- 
ial dark matter halos driven from physical principles. Adopting the picture of 
triaxial collapse based on the Zel'dovich approximation, we derive analytically 
both the minor-to-major and the conditional intermediate-to-major axis ratio 
distributions, and examine how they depend on the halo mass, redshift, and cos- 
mology. Our analytic model is tested against the simulation data given by Jing 
& Suto in 2002, and found to reproduce the conditional intermediate-to- major 
axis ratio distribution successfully and the minor-to-major axis ratio distribution 
approximately. However, the trends of our analytic axis-ratio distributions with 
mass and redshift are opposite to what is found in N-body simulations. This 
failure of our analytic model puts a limitation on analytic approaches based on 
the Lagrangian theory to the halo ellipticity. Given the overall agreement with 
the simulation results, our model provides a new theoretical step toward using 
the axis-ratio distribution of dark halos as a cosmological probe. We also discuss 
several possibilities to improve the model. 

Subject headings: cosmology: theory — large-scale structure of universe galax- 
ies: halos — dark matter 

1. INTRODUCTION 

While shapes of dark matter halos have been conventionally modeled as spherical 
(Navaro, Frenk, & White 1997; Moore et al. 1999), optical, X-ray and lensing observations 
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of galaxy clusters suggest that the shapes of dark halos are far from spherical (West 1989; 
Plionis, Barrow, & Frenck 1991). Recent high-resolution simulations do indicate that the 
density profiles of dark halos are better approximated as triaxial (Frenk et al. 1988; Warren 
et al. 1992; Jing & Suto 2002; Suwa et al. 2003; Kasun & Evrard 2004; Hopkins, Bahcall, & 
Bode 2004). 

Jing & Suto (2002, JS02 hereafter) were able to construct a detailed empirical model 
for the triaxial halo for the first time. Their fitting model turned out to be quite useful in 
quantifying the triaxiality effect on many important observables such as the strong lensing 
arc statistics (Oguri, Lee, & Suto 2003; Oguri & Keeton 2004; Dalai, Holder, & Hennawi 
2004), the halo mass-temperature relation (Yang, Yu, & Shen 2004), and the halo figure 
rotation (Bailin & Steinmetz 2004). 

The most fundamental statistics characterizing the halo triaxiality is the axis-ratio, 
or equivalently, ellipticity distribution functions. Hopkins, Bahcall, & Bode (2004) have 
shown from their N-body simulations that the axis-ratio distributions depend on the halo 
environments in addition to the underlying cosmology. In order to exploit the halo axis-ratio 
distribution as a cosmological probe, therefore, one needs an analytic model beyond the 
empirical fitting formulae. This is exactly what we attempt to propose in this paper. 

In fact, the triaxial shape of a dark halo is a generic prediction of the CDM (Cold Dark 
Matter) paradigm. Bardeen et al. (1986) already derived an analytic expression of the halo 
axis-ratio distribution assuming that dark halos form at peaks of the primordial Gaussian 
density field. They showed theoretically that the CDM dark halos cannot be spherical as 
pointed out earlier by Doroshkevich (1970). However, recent numerical results mentioned 
above have demonstrated that simulated dark halos are even more elliptical than expected 
in the previous density peak approach. This is why we revisit this problem using different 
analytical approaches. 

To construct a new theoretical model for the halo axis-ratio distribution, we adopt the 
cosmic web picture (Bond, Kofman, & Pogosyan 1996) which describes one dimensional 
filaments in the CDM framework. According to the picture, the distribution and spatial 
coherence of initial tidal fields induces the filamentary pattern of the large scale structure. 
Later the filaments bridge between dark matter halos, and the merging of dark halos occur 
preferentially along the bridging filaments. Thus the resulting halos cannot be spherical 
but naturally become elongated along the filaments. The halo ellipticities are expected to 
increase as the hierarchical merging along the filaments proceeds. Therefore, the halo axis- 
ratio can be inferred statistically by combining an evolution model of non-spherical density 
perturbations and the primordial filamentarity of the initial density field. 
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In order to examine the validity of our analytical model, we compare the predictions 
with the simulation results by JS02. For that purpose, we consider two specific sets of 
cosmological parameters that they adopted. The first model is ACDM model which assumes 
that Q m = 0.3, Q\ = 0.7, a 8 = 0.9, and T = 0.2, where Q m and denote the matter density 
parameter and the dimensionless cosmological constant, and <jg is the amplitude of the mass 
fluctuation at 8/i _1 Mpc. Since JS02 use the CDM transfer function of Bardeen et al. (1986) 
neglecting the baryon contribution, the shape parameter T simply equals Vt m h. The second 
model is SCDM model where Vl rn — 1, Q A = 0, a s — 0.55, and T = 0.5. JS02 considered 
dark halos consisting of more than 10 4 particles. Since a mass of a simulation particle is 
2.07 x 10 9 f2 m /i _1 M (512 3 particles in 100/i _1 Mpc comoving cube), it is convenient to define 
the dimensionless mass of M 4 = M/(2.07 x lO 13 fi m /i _1 M ). 

The rest of this paper is organized as follows. Section 2 describes our basic assumptions 
in the analytic modeling. We lay out mathematical details of derivation of the axis ratio 
distributions in §3, and compare the analytic results with their simulation results in §4. 
Finally §5 is devoted to discussion and conclusions. 



2. BASIC ASSUMPTIONS 

To derive the axis-ratio distribution of dark matter halos, we adopt three major as- 
sumptions. First, the trajectory of a dark matter particle in the comoving coordinate is well 
approximated by the Zel'dovich formula (Zeldovich 1970). According to the approximation, 
the key quantities are the three eigenvalues (Ai > A 2 > A3) of the deformation tensor (or 
the tidal shear tensor), dij, which is defined as the second derivative of the perturbation 
potential, at the initial epoch zf 

dij EE didj^l. (1) 

The mapping from the Lagrangian to the Eulerian spaces yields an expression of the particle 
density p in the Eulerian space x at redshift z: 

p(x, z) = = P M = , (2) 

[l- J D + (z)A 1 ][l-D + (z)A 2 ][l- J D + (z)A 3 ]' 

where p is the background density of the universe, and D+(z) = D + (z) / D + (zi) denotes the 
linear growth rate of density fluctuations up to z but normalized to unity at redshift Z{. Note 
that Ai, A2, and A3 in equation (2) represent the eigenvalues of defined at Z{. In practice, 
we use the following fitting formula (Peacock 1999): 

5 r°° 1 + z' 

« + W = 50.10.(1 + *)• + «*]"'/ ly^W^ (3) 
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The linear density contrast, 5, at Zi is simply written as 

<5(x, zj) = Ai + A 2 + A3. (4) 

Second, we assume that a dark matter halo of mass M forms at redshift z when the 
corresponding Lagrangian region (at zi) in the linear density field smoothed over M satisfies 
the following conditions: 

5(zi) = S c (z), X 3 (zi) > X c (z), (5) 

where 5 and A3 are the linear density contrast and the smallest eigenvalue of dij of the 
smoothed density field, respectively. Here 5 c (z) and \ c (z) are redshift-dependent threshold 
of 5 and lower limit of A 3 , respectively: 

5 c (z) = 5 c0 D + (0)/D + (z), X c (z) = X c0 D + (0)/D + (z). (6) 

We use the collapse threshold 5 c q computed in the spherical model. For SCDM, it is 5 C ~ 
1.686. For ACDM, we use the formula given in Appendix B of Kitayama & Suto (1996), 
which depends weakly on cosmology. 

In the spherical approximation, the condition of S(zi) = S c (z) is sufficient for the grav- 
itational collapse at z. In the non-spherical model based on the Zel'dovich approximation, 
however, all regions satisfying 5 = 8 C do not necessarily collapse into dark halos, since A3 can 
be negative even when 5 = 5 C . This is why we impose the additional condition (5) in our 
model based on the Zel'dovich approximation. Nevertheless no reliable modeling is known 
which determines the value of A c0 from physical principles. While the Zel'dovich approxi- 
mation suggests A c0 = 1, those objects are subject to the first-shell crossing, beyond which 
the Zel'dovich approximation is not valid at all. Indeed, Lee & Shandarin (1998) empirically 
proposed a value of A c o = 0.37 in their mass function theory. They argue that realistic 
collapse should proceed along all the three axes almost simultaneously. Thus the collapse 
along the major axis should be accelerated by the collapse along the other two directions, 
resulting in a lower value of A c o. Throughout this paper, we adopt their value of A c o = 0.37 
unless otherwise stated (see §3). 

Third, the principal axes of the inertia tensor of a dark matter halo are aligned with 
that of the linear tidal shear tensor of the corresponding Lagrangian region. Approximating 
that the density profile of a dark matter halo as a triaxial ellipsoid with three distinct axes, 
a, b, and c (we define a < b < c), one can say that the inertia shape tensor of a dark halo has 
three distinct eigenvalues, a,b,c. The three eigenvalues of the halo inertia tensor, {a,b,c}, 
are related to the eigenvalues of the tidal shear tensor, {Ai, A 2 , A 3 } as 

a oc a/1 - D+X 1 , b oc J 1 - D+X 2 , c oc J 1 - D + X 3 , (7) 
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It may be interesting to compare our definition of the halo axes (eq. [7]) with that of the 
density peak formalism (Bardeen et al. 1986). In the density peak formalism, the three 
eigenvalues of the halo inertia tensor are defined as 

111 
aoc— =, bcc—=, coc ^> (8) 

VU VS2 V S3 

where Ci, C2, and (3 are the three eigenvalues of the second derivative of the linear density 
field, didjS. The comparison of two equations (7) and (8) shows that in the density peak 
formalism the inertia shape tensor of a dark halo is almost independent of the tidal shear 
tensor, dij, while in our formalism which basically assumes that the ellipticity of a dark halo 
is induced by the filamentary cosmic web, it is directly related with d^. In fact, the strong 
correlation between the halo inertia and the tidal shear tensors was demonstrated by recent 
N-body simulations (Lee & Pen 2000; Porciani, Dekel, & Hoffman 2002). 

The above three assumptions imply that dark matter halos preferentially form at the 
over-dense nodes of the filamentary web of the initial density field where the principal axes 
of the inertia and the tidal tensors are aligned with each other. Using these simplified 
assumptions, we derive analytically the distribution of the axis-ratios of dark halos in the 
following two sections. 



3. THE MINOR-TO-MAJOR AXIS RATIO DISTRIBUTION 



We start from the joint probability distribution of the three eigenvalues of the tidal 
shear tensor in the primordial Gaussian density field (Doroshkevich 1970): 

3375 / 3/ 2 15/ \ 

p(Ai, A 2 , A 3 ;<7 M ) = —7= — ;-exp ( ^ + ^-)(Xi - A 2 )(A 2 - A 3 )(Ai - A 3 ), (9) 

8V57ro-^ V a M la M J 

where 

h = Ai + A 2 + A 3 , h = AiA 2 + A 2 A 3 + A3A1, (10) 

We define <7 M as the rms fluctuation of the linear density field at Z; L smoothed on mass scale 
M: 

(T 2 M (zi) = J P(k,z t )W* H (kR M )d 3 k, (11) 

where P(k, zi) is the linear power spectrum of the density field, and Wth^Rm) is the top-hat 
filter with R M = [3M/(4vrp)] 1 / 3 . 

We change the variables from {Ai, A 2 , A 3 } to {Ai, A 2 , 5} using equation (4), and find the 
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joint probability distribution of Ai, A2, and 5 from equation (9): 



p(Ai, X 2 ,S;a M ) = 



3375 



8^ 



exp 



7TO" 



M 



35 2 
— r + 



15<J(Ai + A 2 ) 15(A? + AiA 2 + A 



2< 



x(2Ai + A 2 - (J)(Ai - A 2 )(Ai + 2A 2 - 5). 
Applying the Bayes theorem and the Gaussian distribution of the linear density 



p(S; <Jm) = 



exp 



'hi 



the conditional probability distribution of Ai and A 2 at 5 = 5 C is written as 

p(\ 1 ,X 2 ,S = 5 c ;a M )dS 



p(Ai, X 2 \5 = 5 c ;a M ) = 



p(S = 5 C ; a M )d5 



(12) 



(13) 



exp 



55 2 c 155 c (Ai + A 2 ) 15(A? + AiA 2 + A|) 



3375^ 

x(2Ai + A 2 - (J c )(Ai - A 2 )(Ai + 2A 2 - 5 C 



2< 



2< 



2< 



(14) 



We define two real variables, /ii and /i 2 , as the axis ratios of a triaxial halo: 

b a 



/ii = 



A*2 = 



(15) 



(/i 2 < /ii < 1). According to the third assumption in §2, /ii and /i 2 are written in terms of 
Ai, A 2 , and 5 as 



1 - D . A 5 



'' ! " 1 A j - A 2 ) 
with the following constraints: 



1^2 



1 - D + Xi 



1 - £>+(5 c - Ai - A 2 



(16) 



Ai - A 2 > A c (z), Ai < -~ 



1 



L> + (z) 



(17) 



The first constraint A3 > \ c (z) guarantees the collapse along all three axes, in accordance 
with the second assumption in §2. The second constraint Ai < l/D + (z) guarantees that 
Hi and H2 are all real. Note that if Ai < l/D + (z), then it automatically implies A3 < 
A 2 < l/D + (z). According to equation (14), however, Ai has a non-zero probability of 
Ai > l/D + (z). We simply do not consider the parameter region of Ai > l/D + (z) since they 
correspond to the break-down of the Zel'dovich approximation in the non-linear regime after 
the first-shell crossing (Ai = l/D + (z)). 



-7- 



Now, we write the probability density that a dark matter halo of mass M formed at 
redshift Zf has a intermediate-to-major axis ratio of b/c and the minor-to-major axis ratio 
of a/c as 

p(b/c,a/c;M,z f ) = fj, 2 \S = S c ;a M ; zf) 



= Ap(X 1 , X 2 \8 = 5 c ;a M ;z f )Q\ Ai 



x0[5 c -A c -(A 1 + A 2 ) 

where we solve equation (16) for Ai, A2 and A3 as 

l + (D f 6 c -2)i4 + i4 



(d\ 1 d\ 2 ) 



(dfiidfi 2 ) 



(18) 



Ai = 
A 2 = 
A3 = 



l + (D f 5 c -2) f i 2 1 + f 4 
D f 5 c -2 + yl + 



(19) 
(20) 
(21) 



In the above expression, 5 C and Df depend on the formation redshift, zf. S c (zf) and Df 
D + (zf), and 6 is the Heaviside step function. The normalization constant A satisfies 



J p(Xi, X 2 \5 = S c ; (7 M , z f )9 - A^ Q[5 C - A c - (Ai + Xfjjd^d^ = I. 



(22) 



Finally we denote by \{dXidX 2 )/{d^Lid^ 2 )\ the Jacobian of the transformation from {Ai, A 2 } 
to {fii,fi 2 }, which we find from equations (19) and (20): 



(<9Ak9A 2 



(<9/ii<9/i 2 ) 



4(D f 5 c -3) 2 ^ 2 



(23) 



Integrating equation (18) over b/c, we find the probability density that a dark halo of 
mass M formed at redshift Zf has a minor-to-major axis ratio of a/c: 



p(a/c;M,z f ) = I p(b/c, a/c; M, Zf)d(b/c) 

J a/c 



p[fJ>i, /J-2\S = S c (z f ); a M \dn\ 



I Ap[\ 1 ,\ 2 \6 = 6 c (zf);<T M ]Gl-~--\ 1 ] 
h2 \Df J 
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x0[5 c -A c -(A 1 + A 2 )] 



(<9Ai(9A 2 ) 
(djiidji2) 



djii. 



(24) 



Equation (24) is the axis ratio distribution at the formation epoch, Zf. The axis-ratio 
distribution at the observation epoch, z, can be readily found as 



where the formation epoch distribution, dpf/dzf, represents the probability that a halo 
of mass 2M that exists at z had a mass greater than M for the first time at Zf. Since 
the formation epoch distribution in the current non-spherical model is almost impossible to 
work out analytically, we use the spherical counterpart instead. In practice we use the fitting 
formula by Kitayama & Suto (1996) for the analytic expression derived by Lacey & Cole 



Let us emphasize the difference between equations (24) and (25): while the former gives 
the probability density that a dark halo of mass M has a minor-to- major axis ratio a/c at 
its formation redshift Zf, the latter is the counterpart evaluated at the observation redshift 
z(< Zf). Of course, in numerical simulations and observations, the latter, not the former, is 
the relevant observable. In what follows, therefore, we mainly consider the latter evaluated 
at the observation epoch, z. 

Figure 1 plots the A c0 -dependence of p(a/c; M, z) in the ACDM model: A c0 = 0, 0.1, 0.2, 
0.3, 0.35 and 0.4 (dotted, dashed, long-dashed, solid, dot-dashed, and dot-long-dashed lines, 
respectively). For this plot, we choose the halo mass M4 = 10 at z = for definiteness. The 
axis-ratio distribution of dark halos shifts toward the high axis-ratio side (more spherical) as 
A c o increases. This is theoretically understandable since higher values of A c o correspond to 
those Lagrangian regions whose major axis lengths are closer to the other two (see eq.[7]). 
Given such, we adopt A c0 = 0.37 in the analysis below, following Lee & Shandarin (1998) 

Figure 2 shows how p{a/c; M, z) depends on the halo mass and the redshift for the case 
of ACDM. The upper panel shows the z-dependence of p(a/c; M,z) for a fixed mass scale 
of M 4 = 10 at redshift z — 0, 0.5, 1, 1.5 and 2 (solid, dotted, dashed, long-dashed, and 
dot-dashed) respectively. While the lower panel plots the mass-dependence of p(a/c;M,z) 
at z = for M 4 = 1, 5, 10, 20 and 30 (solid, dotted, dashed, long-dashed, and dot-dashed) 
respectively. As can be seen, the analytic distribution p(a/c; M, z) depends on the halo mass 
and redshift consistently: the distribution moves toward the high axis-ratio section as the 
halo mass and the redshift increase. 

In other words, our analytic model predicts that the more massive a dark halos is, the 
less elliptical it is, and that a dark halo of given mass is less elliptical at earlier epochs. These 




(25) 



(1994). 
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two trends are consistent with the theoretical work of Bernardeau (1994). By applying the 
perturbation theory to a primordial Gaussian density field, he proved theoretically that the 
larger halos formed at higher density peaks are rounder. However, these trends are opposite 
to what is found through accurate calculations of N-body simulations: It was found in 
simulations that the more massive halos are more elliptical, and that the halos of given mass 
observed at earlier epochs are more elliptical (Bullock 2002; Springel, White, & Hernquist 
2004; Jing & Suto 2002; Hopkins, Bahcall, & Bode 2004). This failure of our analytic model 
implies that the dependences of the axis-ratio distributions on mass and redshift are not 
fully determined by simply applying the Zel'dovich approximation. 

Figure 3 shows our analytic prediction (eq. [25]) at z — 0, 0.5, and 1 (top, middle, and 
bottom, respectively) for three mass scales M 4 = 1,4, 10 (solid, dotted, and dashed, respec- 
tively). They should be compared with the numerical results (histograms) from JS02 for the 
ACDM (left panels) and the SCDM (right panels). Clearly the analytic predictions agree 
with the numerical results reasonably well; they reproduce well shapes and characteristic be- 
haviors (especially for ACDM) such as the peak positions, the dispersions, and the decrease 
of the mean axis-ratios with the increase of redshift. On the other hand, we notice that 
the numerical histograms slightly move toward the low axis ratio section as the halo mass 
increases, which disagrees with the analytical predictions. We discuss on this disagreement 
between the analytical and the numerical results in §5. 



4. THE CONDITIONAL INTERMEDIATE-TO-MAJOR AXIS RATIO 

DISTRIBUTION 

The probability density, p(b/c; M; z), that a dark halo of mass M is observed at redshift 
z to have an intermediate-to-major axis ratio of b/c can be also computed in a similar 
manner. However, to investigate the overall triaxiality of a dark halo, what is more relevant 
is the conditional probability density distribution, p(b/c\a/c; M; z), that a dark halo of mass 
M is observed at z to have an intermediate-to-major axis ratio b/c provided it has a minor- 
to-major axis ratio a/c. In principle, one can find this conditional probability density from 
the Bayes theorem: 

/ 1 i j , , p(b/c,a/c: M: z)d(a/c) . . 

p(b/c\a c;M;z) = FK ' -. 1 ' \ ' \ ; . 26 
FK ' ' 1 ' p(a/c;M;z)d(a/c) y ' 

Although it is possible to derive p(b/c\a/c; M; z) analytically as well, it is not easy to con- 
struct the statistical sample either from the current simulations or from observations; there 
would be only a few dark halos of mass M at redshift z with the fixed minor-to- major a/c. 



To overcome the poor number statistics, JS02 combined all the halo mass, i.e., they 
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computed p(b/c\a/c; z) instead of p(b/c\a/c; M; z). For a direct comparison with their result, 
we compute p(b/c\a/c; z) according to 

p(b/c\a/c; z) = I dM f°° dz f n{2M; z ) dpf ^ 2M ' Z l p(b/c\a/c; 2M; z f ), (27) 

JSm Jz OZf 

where p(b/c\a/c;2M; zf) is given as equation (18), n(M;z) represents the number density 
of dark halos of mass M that exist at z, and Sm represents the whole mass range to be 
considered. In practice, we use the fitting formula given by Sheth & Tormen (1999) for 
n(M; z). 

Figure 4 illustrates the redshift-dependence of p(b/c\a/c; z) in the ACDM model at 
z — 0, 0.5, 1, 1.5 and 2 (solid, dotted, dashed, long-dashed, and dot-dashed) respectively. 
The minor-to-major axis ratio is fixed to be a/c = 0.55. As one can see, the conditional 
distribution p(b/c\a/c;z) is insensitive to the redshift, which is also consistent with the 
finding of JS02. 

Figure 5 compares the analytic predictions (curves) with the numerical findings (his- 
tograms) for a/c = 0.55, 0.65, and 0.75 (top, middle, and bottom, respectively) at z = 
in ACDM (left panels) and SCDM (right panels) models. The histograms are computed 
by averaging over 0.5 < a/c < 0.6, 0.6 < a/c < 0.7, and 0.7 < a/c < 0.8 (top, middle, 
and bottom, respectively). The agreement between analytic and numerical results is quite 
satisfactory. 

5. DISCUSSIONS AND CONCLUSIONS 

We have derived an analytic expression for the axis-ratio distribution of triaxial dark 
matter halos for the first time. In constructing the analytic model, we adopted the cosmic 
web picture in which the ellipticities of dark matter halos are induced by the coherent tidal 
fields in the initial density fluctuations, and employed the Zel'dovich-type collapse condition 
as a diagnostics. Our analytic model is successful in reproducing the basic behaviors on a 
qualitative level found from the previous numerical simulations (JS02). This ensures our 
basic picture that the origin of the halo triaxiality is the filamentarity of the initial density 
field. Nevertheless we found a discrepancy with the numerical results on a quantitative level; 
in particular the predicted dependences of the halo triaxiality on mass and redshift seem 
inconsistent with the numerical results. 

In order for the future improvements, let us critically discuss possible caveats in our 
analytic approaches: First of all, we use a cooked-up collapse condition of 5 = 5 C and 
A3 > A c . What we really mean to have is a practical and simple collapse condition for 
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the formation of a triaxial dark halo in the nodes of filamentary web of the density field 
by combining the density peak formalism and the Zel'dovich approximation. This collapse 
condition is theoretically unjustified. 

How to identify a triaxial halo and where to locate the corresponding Lagrangian site 
in the linear density field is a touchy issue in the non-spherical dynamical model. Although 
Bond & Myers (1996) proposed the peak patch picture as a complete non-spherical model, 
their formalism is too complicated to follow analytically in practice. Besides, as discussed 
in Suto (2002), there is still no unanimous agreement among theory, observations, and sim- 
ulations about how to define a gravitationally bound object in practice. The disagreement 
among the three become more serious when a dark halo is to be described as triaxial. Hence, 
to find a theoretically justified criterion for the formation of a triaxial dark halo and to derive 
the axis-ratio distribution more rigorously with the criterion, it will be necessary to address 
this difficult issue first, which is beyond the scope of this paper. 

As Bernardeau (1994) proved theoretically that a rare event in the linear density field 
is inclined to be quite spherical. In other words, a linear over-dense region of high-mass 
must be more spherical than a low-mass over-dense region. This is the case that our analytic 
distribution predicts. The difference in the trends with mass and redshift between our 
model and the simulation results may be caused by the fact that in reality the shape of a 
dark halo must be affected by subsequent nonlinear clustering process. We argue here that 
the nonlinear merging event must play a key role in increasing the halo ellipticity. Many 
N-body studies demonstrated that the merging event occurs anisotropically along filaments 
(e.g., West, Villumsen, & Dekel 1991; van Haarlem & van de Weygaert 1993; Dubinski 1998; 
Faltenbacher et al. 2002). This anisotropic merging event tends to make the shapes of dark 
halos elongated by aligning their substructures along with the orientation of their major 
axes (e.g., Knebe et al. 2004). Moreover, it was also shown that the merging process is more 
rapid for the case of higher mass halos (Zhao et al. 2003). The tendency of being more 
spherical in the higher-mass section of the linear density field (Bernardeau 1994) is likely to 
be compensated by the anisotropic merging effect. Therefore, the overall dependence of the 
axis-ratio distribution on the halo mass and redshift may be weakened by this compensating 
effect. Our future work is in the direction of incorporating semi- analytically the anisotropic 
merging process into our analytic model. 

In passing, it is interesting to note that the anisotropic merging of dark halos along 
filaments has an implication about the mass function of dark halos. In the standard mass 
function theory based on the Press-Schechter theory, the mass function is independent of 
the power-spectrum. If the merging really occurs in an anisotropic way along filaments, 
however, the filamentarity in the medium must affect the final mass distribution of dark 
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halos. For example, for the case of a power-law spectrum, the mass function might depend 
on the power-law index sensitively, since the filamentarity of the density field should depend 
on the power-law index. We wish to present the effect of anisotropic merging on the mass 
distribution of dark halos in the near future (Lee & Jing 2005, in preparation). 

The discrepancy between our analytic model and the simulation results on the mass and 
redshift dependences implies that the shapes of dark halos are not fully determined by simply 
applying the Zel'dovich approximation. Rather, one has to take into account complicated 
nonlinear gravitational clustering like the anisotropic merging event. Nevertheless, the overall 
agreement of our analytic model with the simulation results gives us a hope that our analytic 
model will be useful in quantifying how the initial cosmic web induces the ellipticity of dark 
matter halos, and provide a first theoretical step toward the goal of using the ellipticity 
distribution of galaxy clusters as a new cosmological probe (Lee 2005, in preparation). 
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Fig. 1. — Probability density distribution of the axis ratio a/c of halos for the five different 
cases of the short-axis cut-off in the A CDM model; A c0 = 0, 0.1, 0.2, 0.3, 0.35 and 0.4 
(dotted, dashed, long-dashed, solid, dot-dashed, and dot-long dashed lines) respectively. 
Here, the halo mass and redshift are set to be M 4 = M/(2.07 x lO 13 O m /i -1 M ) = 10 and 
z = 0, respectively. 
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Fig. 2. — Behavior of the probability density distributions of the axis ratio a/c of halos 
with the change of observation epochs Zf and the halo mass M in the A CDM model. Top: 
Case in which z = 0, 0.5, 1, 1.5 and 2 (solid, dotted, dashed, long-dashed, and dot-dashed, 
respectively). Bottom: Case in which M 4 = 1, 5, 10, 20 and 30 (solid, dotted, dashed, 
long-dashed, and dot-dashed,respectively), where M 4 = M/(2.07 x lO 13 f2 m /i _1 M ). 
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Fig. 3. — Probability density distributions of the axis ratio a/c of halos in the A CDM 
(left) and SCDM (right) models at three different observation epochs; z = 0, 0.5, and 1 (top, 
middle, and bottom panels, respectively), with the choice of A c0 = 0.37. In each panel, 
the histograms represent the numerical results (Jing & Suto 2002) for those halos of mass; 
1 < M 4 < 2, 2 < M 4 < 6, and M 4 > 6 (thin solid, dotted, and dashed lines, respectively.) 
where M 4 = M/(2.07 x lO 13 f2 m /i _1 M0), while the curves represent the analytic results (see 
eq.[24]) for those halos of mass; M 4 = 1, M 4 = 4, and M 4 = 10 (thick solid, dotted, and 
dashed lines), respectively. 
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Fig. 4. — Behavior of the conditional probability density distributions of the axis ratio b/c 
at z — 0, 0.5, 1, 1.5 and 2 (solid, dotted, dashed, long-dashed, and dot-dashed), respectively. 
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Fig. 5. — Conditional probability density distributions of the axis ratio a/c of halos in the 
A CDM (left) and SCDM (right) models provided that the minor-to- major axis ratio has a 
certain value: a/c = 0.55, 0.65 and a/c = 0.75 (top, middle, and bottom panels), respectively 
at z = 0. In each panel, the histogram and the curve represent the numerical (Jing & Suto 
2002) and the analytic results, respectively. 



